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The aim of this work is to solve the dispersion relations near the first excitation threshold of photon propagating along the 
magnetic field in the strong field limit. We have calculated the time damping of the photon in two particular cases: the 
degenerate gas as well as the diluted gas limit being both important from the Astrophysical point of view. In particular 
the diluted gas limit could describe the magnethosphere of neutron stars. The solutions have been used to obtain a finite 
Quantum Earaday angle in both limits. A resonant behavior for the Earaday angle is also obtained. To reproduce the 
semi-classical result for the Earaday rotation angle the weak field limit is considered. 
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1 Introduction 

Photon propagating in a magnetized medium gives rise to 
the existence of different dispersion relations which lead to 
many observable effects. In the Astrophysical and Cosmo¬ 
logical scenarios one of the most important ef fects is the 
Faraday rotation (FR) angle ( Giovannhiil 1997h which al¬ 
lows to estimate values of the magnetic field in the Universe. 

The issue of the light propagation in a material medium 
like insulators, metals etc has been covered extensively 
from both the classical and semi-classical point of view. 
However this approach is not valid anymore when the 
light travels from inside a compact object characterized 
by huge densities and strong magnetic field. An approach 
from the perspective of Quantum relativistic theor y shoul d 
be included in these stud ies ( Rojas & ShabadI ( 1982 )- 
Cruz Rodriguez et al.l ( 120 13b and references therein). 

In a previous paper (ICruz Rodriguez et al.l (12013^ 1 the 
close relation between Quantum Faraday Effect (FE) and 
Quantum Hall Effect was tackled. The starting point of that 
study was the solution of the Maxwell equation which re¬ 
sults from considering the radiative correction: magnetized 
dense self-energy of photons. The photon self-energy gives 
count of the propagation of the photon through the magne¬ 
tized medium (electron-positron plasma). 

Nevertheless there are no observations capable of mea¬ 
suring neither the Earaday angle nor the value of the mag- 
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netic field from light that leaves the inner regions of the 
compact stars. There is also no terrestrial experiments de¬ 
signed to study a similar scenario yet. However, there are 
experimen ts where the ER a ngle has been measured in 
Graphene ( Crassee et al.l201l|) , so it would serve as a test of 


our res ults. Hence, our previous paner lCruz Rodriguez et al 


( 2013b had also the purpose to obtain the quantum Earaday 
angle and QHE in 2 Dh- 1 systems (Graphene) from the 3 Dh- 1 
results. 

The aim of the present paper is to continue the stud- 
ies of the propagation of a photon in magnetized medium 
( Cruz Rodriguez et'H] 201 31 2014 ) thinking in the Astro- 
physical implications. In particular we devote this work to 
solve the dispersion relation near the first excitation thresh¬ 
old. The behavior near the thresholds are related to absorp¬ 
tion processes and a complete description of the dispersion 
law of photon requieres to consider them. For instance, in 
our first calculation of the Quantum Faraday angle (for a de- 
generate electron-positron p lasma at strong magnetic field) 
( Cruz Rodriguez et'^l2013b we obtained a non-finite value 
of the angle due to the branching points associated with the 
absorption processes. As we will see below these singular¬ 
ities can be removed finding the solutions near the singular 
points. Qur study is done in two physical scenarios which 
are interesting for Astrophysics: degenerate gas (/i ^ T) 
and diluted gas {p, < T). For the diluted gas limit we con¬ 
sider both: strong (eB ^ p?) and weak (eB <C p^) field 
regime. 

The paper is divided as follows. In Sec |2]the solution 
of the dispersion equation near the first excitation threshold 
is solved and the Faraday angle obtained in the degenerate 
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and strong field limit is corrected considering the absorp¬ 
tion. Section [3 is devoted to compute the scalars r and t in 
the diluted gas approximation in order to solve the disper¬ 
sion equation and to hnd the Faraday angle. The weak field 
limit is also studied in the diluted gas approximation with 
the aim of reproducing the semiclassical results. Finally, in 
Secl4]we state the concluding remarks. 


2 Solution of the dispersion equation near 
the first excitation threshold 


The propagation of photons in a relativistic medium at h- 
nite temperature and density, in the presence of a constan t 
magnetic field was studied in ( Roias & Shabadll9^ , 1982 ). 
The Maxwell equations in the Fourier space are 

[{k'^g^Lu - k^K + = 0, (1) 

where a^, are small corrections to the photon field ampli¬ 
tude, zl® is the external magnetite held and 11^;^ is the pho¬ 
ton self energy which contains all the information related to 
the interaction with the medium {v, fjL = 1, 2, 3,4). 

With the aim of solving the dispersion equation in order 
to hnd the photon propagation modes, the self energy ten¬ 
sor is di agonalized and a dispe rsion law for each mode is 
derived ( Roias & Shabad ( 1982h f. 


= Ki{k), 


( 2 ) 


where, /cp -F-F /c^=A:| -F A:| = k'^ — w'^. 

For the particular case of propagation along the mag¬ 
netic held {k± = 0), in the charge asymmetric case, there 
are three non-vanishing eigenvalues, one is a pure longitudi¬ 
nal mode and the other two are transverse modes with differ¬ 
ent di spersion laws, which is the key to the Faraday Effect 
(FE) dRoias &Shabadl(ll982h '). 

The eigenvalues related to the transversal modes could 
be written as k,i ,2 = ir = t ± Ir, which give rise to the 
dispersion equation 


where. 


kfi = kl—u}"^ = ki,2, 


r = il,, t = -^Iu 


(3) 


and Ir, It are the integrals 

e^Buj 


I. =- 


^;^/t^F3/(F3, ks, B,u},fj,){neiep^n) - np{sp^n)), (4) 


( 3 ) / + 2eB(n -F n') 

/(p3,fc||,/i,S,T) = ^Bf2,(0) J - - - 


n,n' 


D = [2psk3 + A:j| + 2eB{n - n)f - 

oo 

It — J' /b||,/r, B, T)(?T.e(£rp yj)-F (5) 


^=E 

n^n' 




(2p3A:3 -F Jnn' )(fcy -F 2eB(n -F n')) 
D 


where ne,p(ep,n)=lXl + are the Eermi-Dirac 

distribution for electrons and positrons, n is the Landau 
quantum number, the energy levels are given by „ = 
-F -F 2neB and: 

Jnn>=k‘\\+2eB{n'-n), 

Due to the singular denominators in (|4|i and (|5]l, the in¬ 
tegral can be divided in a real and imaginary contribution 
(Perez & Shabad 1982), 


1 

s — u! + ie 


PV- 


S — (jJ 


iTreS{s — cu). 


( 6 ) 


Ear from the threshold only the contribution of the prin¬ 
cipal value of the integral is relevant, but near the singulari¬ 
ties the imaginary contribution should be considered. 

In Cruz Rodriguez et"^ ( 2013 ) an approximate solution 
of the dispersion equation was derived, regarding the depen¬ 
dence of the scalars r and t with k^ (long wave limit). In par¬ 
ticular the strong held (eB>>/i^) and degenerate {^^>>T) 
limits were considered due to its relevance in astrophysics. 
However, only the contribution of the PV of the integral was 
included. Nevertheless, a complete study of the dispersion 
equation should include the solution near the threshold. 

Close to the branching points the dispersion equation 
takes the form: 


kfi = Im{t) ± Im{Ir). 


(7) 


The hrst excitation threshold which give rise to an inhnite 
behavior is n = 1 and n = 0 or vice versa. Around this 
branching point t he dis persion relation can be written as 
( Roias & Shabad ( 1982h l: 


kt, = 


e-^B 


2eB 


Att (fcy — -F 2eB)i/'^ 


n: 


( 8 ) 

where (fc^ = — [m — (m^ -F 2eB)^/^]^). In the degenerate 
limit (/3 —>■ c») the distribution is the difference between 
two step functions: 

NeiSq) = - Sq) “ - Eg -Uj). 

After squaring (|8]l we may write it as a cubic equation: 


+ px^ + g — 0, 

where x = kfJrn^, is a non-dimensional variable and 


(9) 


kf 


p=-- 


q=- 


m 


'e^B(kf -F2eB)' 


4(to^ -F 2eB)i/^ 


dTrm"^ 


( 10 ) 
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Upon the variable substitution y = x — | the square term is 
removed ( Birkhoff & Mac Lanai 1965 ). 


y'^ + Ay + B = {), 

A= ^(3q-p^), B = ^(2/-9m + 27), 

and the discriminant is defined as: A = 127 + 5^/4, if 

A is greater than zero the solutions are: 


yi=M + 7V, y2.3 = -i(M + iV)±i^(M-iV), 

( 11 ) 


Above the threshold the imaginary part of the solution is 
the responsible for the time damping of the photon e*^*^*! 

gi(i^+ir)t Qa^sed by the excitation of electrons to higher 
quantum states. 

From the equation w+iF = (A:| —x| we obtain 

the imaginary contribution to the photon energy (Hugo et al. 
1982): 


r =—|(m^i?e[x2,3])^ +("i^/to[x 2,3])^|= -|w^i?e[x2,3] |^ ■ 
v2 

( 12 ) 

r is the probability of absorption of the eingenmodes which 
remains finite at the threshold and depend only on the value 
of the magnetic field. In the next section we are going to use 
this contribut ion in order to continue the p revious study of 
the FR angle ( Cruz Rodriguez et al. ( 20131) ). 


2.1 Faraday Effect in the degenerate limit 


In the study initiated in ICruz Rodriguez et al.l (120131) the FR 
angle, 

(fc_ — k+)L 


e = 


(13) 


was obtained from the solutions of the dispersion equation 
(O in the long wave limit (ks 0), regarding the contri¬ 
bution near the branching point. But at this point, using the 
result derived in the above section, F can be added to the to¬ 
tal photon energy in the expression obtained for the Faraday 
angle. 


e ^ ir{nuj,B) 9 ^ ir{nu} + ir,B) 

L 2hu!c L 2h{ijj-\-iT)c 


(14) 


This correction was consider ed for graphene-like systems in 


Cruz Rodriguez et al.l (12013h , including it as a phenomeno¬ 


logical parameter. But now we are able to estimate the con¬ 
tribution of the absorption process using (fT2l) for our 3 Dh- 1 
fermion gas. 

As can be seen in FiglT] a finite angle is obtained. The 
enhancement of the Faraday angle is associated with the cy¬ 
clotron resonance, the curve turn smoothed due to the con¬ 
tribution of the absorption term. The inset figure shows the 



Fig. 1 (Color on line) FR angle as a function of the photon 
energy is plotted. Typical values of the magnetic field and 
photon energy are considered: Tiuj from IMeV to fiOMeU, 
B — The value estimated for F is 13MeV. 


dependence on B for IMev and 20MeV, which is approx¬ 
imately lineal. 

This result could be compare to the experimental mea- 
surements of FR an gle in graphene grown over a substrate 
( Crassee et alj201ll) which behaves as a tridimensional sys¬ 
tem. As we can see Fig[T]shows a go od agreement with th e 
behavior of the FR angle reported bv ICrassee et al.l (1201 ih . 


3 Dispersion equation in the diluted gas limit. 


Another interesting limit from the astrophysical point of 
view is the photon propagation in a relativistic and diluted 
media (T >> /i), where strong magnetic field exist. This 
is the case of radiation propagating through the magnetho- 
sphere of neutron stars. 

If (T > > /i) the particles distributes classically. 


_ 1 
“ 1 _|_ g(EeT/i)/3 




(15) 


After a Taylor expansion of the functions / and p in (01 and 
(|5]), using the fact that the main contribution to the integrals 
arise from the P 3 = 0 term, so the integrals Ij. and It can be 
written as: 


Ir 


eBoj 

2tt'^ 


f{0,k3,uj,y,B,T) 


dpsiUe - Up), 


(16) 


It = g{0,k3,uj,p,B,T) j dp3{ne + np). (17) 

By expressing the integral in terms of the modified Bessel 
function, we get 

Arrr 

dp3ne,p = - ^ —A:i(to„/3)F(3/2), (18) 

where 



K^{z) 


yfHz'' 

2''T{v+l/2) 





ly-^/^dt, (19) 
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and TTin = \/m? + 2neB. 

In the strong field limit only the contribution of the low¬ 
est Landau level (LLL) remains, which means to consider 
in (IHl and (|5]l, n = 0, n = 1 or vice versa, depending on 
the Kronecker delta expressions. Finally in the diluted and 
strong field limit, 

8e^B + 2eB)u}mKi{ml3)T{3/2) sinh(/i/3) 
7r®/2 (fcy -I- 2eBY — 

( 20 ) 


8e^B 2w'^rin?Ki{ml3)T{i/2) cosh(/r/3) 
TT®/^ (/cjl -I- 2eBY — 


and from OJ, the solution of the dispersion equation is 
found, 

0^3 D 

k] = ^iCi(m/3)r(3/2) 


2uj^m^ cosh(^/3) ± muj(kn + 2eB) sinh(/r/3) 


y (fcy -I- 2eB)'^ — 

Holding in mind the condition (T >> /r). 


( 22 ) 


g/i/? 

sinh(/i/3) ~ cosh(/r/3) ~ 

and as the Bessel function in (1221) converges quickly to zero, 
the dispersion equation becomes quadratic for each mode, 
which solution takes the form; 


=---h \/{2eB ± 2mui)‘^ — 2A{B, T, fj,)muj, 

(23) 

where A(B,r,/r) = ^iFi(m/ 3 )r( 3 / 2 )e^' 5 _ 

Then, following the same procedure as in the above sec¬ 
tion, the solution near the threshold is studied. Taking as a 
starting point the dispersion equation near the first threshold 
dUl, but with the proper distribution in the diluted approxi¬ 
mation, the value of F is estimated. 

In the next section the Faraday angle in the diluted gas is 
studied considering also the solution near the first excitation 
threshold. 


3.1 Faraday angle in the diluted gas limit. 

From the general definition of the Faraday rotation angle 
(HI, where k- and k+ are the wave vectors of each mode 
obtained from (| 2 ^ . the following result is derived (units are 
recovered); 

^ g{eBtic)tiui 

L 8{{eBhcY - {Tiw fy ^ ’ 

g{eBhc) = 



Fig. 2 (Color on line) FR angle versus photon energy is 
plotted. The photon energy Tiuj runs between O.lMeL to 
lOMeL, the magnetic field B = lO^^G and F = l.SMeV. 


the temperature. The singularities in the denominator are 
solved adding the contribution of absorption calculated in 
the previous section. 

The parameters used to evaluate the Faraday angle and 
the absorption coefficient, correspond to typical values of 
magnetic field in the magnethosphere of neutron stars (B = 
lO^^G) and with radiation emitted in the 7 region of the 
spectrum, where the denominator in (l24l i is close to the sin¬ 
gularity. 


9 g{eB?ic){Ti(jj + iT) 1 

L^~ 8 {eBncy-{nuj + iTy\ ’ 

(25) 

In Fig|2]the Faraday angle is plotted as a function of the 
photon energy in units of the maximum value of 

the angle is close to the region were hui = eBfic and it is 
related to electronic transitions to higher Landau levels due 
to photon absorption. 


3.2 Weak Field Limit 


Starting from (fl4l i our aim now is to study the weak field 
and diluted gas limit (/r^ < w? <C eB <C ) in order to 
reproduce the semi-classical results for the FE. 

In this limit the energy separation between two con¬ 
secutive energy levels (quantized by the Landau levels) 
is proportional to VeS and the spectrum recovers e = 

-I- TOp (p^ = p| -I- p5_, p 5_ = 2|ei?|n) so only the 
n = 0 contributes to the sum. The rest of the sum over Lan¬ 
dau levels higher than cero n y 0 can be replaced by the 
integration over p±, 


\eB\ f (Pp± 

27r J (27r)2 ' 


(26) 


In the diluted limit the Faraday rotation depends on the 
photon energy, the magnetic field, the particle density and 


Upon this substitution and considering the classical 
limit for the particle distribution we have now an integral 
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over the whole momentum space, 


It = ^sinh(/r/3)^y (fp 


J + 


+ m^) 


- / d^P 


J + 2p‘] 


J '2 — 4a;^(p2 + TO^) 


= -e/3 


(27) 


where = p§ + p^.’ + 2ei7, J = —— 2ei3 

and = p^ + In the relativistic limit ~ p^ and 
the integral may be written in terms of the energy (dPp = 
Ap^dp = ‘ie^de) 


Ir = sinh(p/3) 


de 


(J+|eV 

J2 _ 4^2£2 

(/ + |e 


— de 


J'2 - 4a;2e2 ^ ’ 


(28) 


Using the dehnition of the function Fn{z) (D’Olivo et al. 
2003), eq. (1281) can be written as 


Ir = sinh(p/ 3 ) (^J F2{z )- JF2{z)^ 

where z = z' = ^ and 


Fr^iz) 


r°° x"e-“ 



(30) 


Ei{z) is the exponential integral function. 

F 2 {z) = 1 - ^{e-^E,{z) - e^E,{-z)), 

and 

E^iz) = 2 - ^{e-^E,{z) - e^E,i-z)). 

Expanding around z = 0, taken into account only the terms 
of first order in B, using the dehnition for the Faraday angle 
given in (Ell, the classical limit is derived. 


Op ^ 8e^B 


(31) 


where also the approximation sinh(p/3) ~ p/3 was used. 
This is the well known semi-classical result for the Faraday 
rotation angle (Ganguly et al. 1999 & D’Olivo et al. 2003). 


4 Conclusions 

We have improved the study of the dispersion relation of 
photon propagating along the magnetic held. We have ob¬ 
tained the solutions near the hrst excitation threshold. The 
study is general but has been applied in two limit cases: the 
degenerate gas (p ^ T) and diluted gas (p ^ T) due to 
their relevance in Astrophysical scenario. Our outcome are 
the following: 


- The result of the damping time F related to absorption 
process is general (for the hrst excitation threshold) and 
it can be used in any particular limit. 

- FR angle can be obtained as a hnite quantity after in¬ 
troducing r in the singular denominators. On the other 
hand FR angle continues showing a resonant behavior. 

- The behavior of the FR angle with the energy in the de¬ 
generate limit has a direct analogy with 
porte d for Graphene grown over substrate 

201 lh . 

A future perspective of this work would be to apply our 
study to the magnetosphere of the neutron stars described 
by the diluted gas. The scattering and absorbtion processes 
of the radiation in that region of the neutron stars could an¬ 
swer some open question: as cooling, braking index etc of 
neutron stars. Our calculation of the damping time would be 
useful to develop this task. 
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